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Abstract: We introduce a model of a randomly growing interface in multi¬ 
dimensional Euclidean space. The growth model incorporates a random or¬ 
der model as an ingredient of its graphical construction, in a way that repli¬ 
cates the connection between the planar increasing sequences model and the 
one-dimensional Hammersley process. We prove a hydrodynamic limit for the 
height process, and a limit which says that certain perturbations of the ran¬ 
dom surface follow the characteristics of the macroscopic equation. By virtue 
of the space-time Poissonian construction, we know the macroscopic velocity 
function explicitly up to a constant factor. 


1. Introduction 


We introduce a model of a randomly growing interface, whose construction in¬ 
volves the height of a random partial order. The interface is defined by a height 
function on d-dimensional Euclidean space, and the related model of random order 
is in d + 1 dimensional space-time. Our goal is to emulate in higher dimensions 
the fruitful relationship between the one-dimensional Hammersley process and the 
model of increasing sequences among planar Poisson points. The connection be¬ 
tween Hammersley’s process and increasing sequences was suggested in Hammers¬ 
ley’s paper J^| , first utilized by Aldous and Diaconis [lj , and subsequently in papers 
2 1| . |22I . I25l . |27| . A review of the use of Hammersley’s process to study increasing 
sequences appeared in 0 , and of the wider mathematical context in [3] • The study 
of higher dimensional random orders was started by Winkler [30]. 

The interface process we introduce is defined through a graphical representation 
which utilizes a homogeneous space-time Poisson point process, and in particular 
the heights of the partial orders among the Poisson points in space-time rectangles. 
This definition suggests a natural infinitesimal description, which we verify in a 
sense. After defining the process, we prove a hydrodynamic limit for the height 
function. This proceeds in a familiar way, by the path level variational formulation. 
The deterministic limiting height is the solution of a Hamilton-Jacobi equation 
given by a Hopf-Lax formula. 

Next we use this process to prove a limit that in a way generalizes the law of 
large numbers of a second class particle in one-dimensional systems. In interacting 
particle systems, a second class particle is the location X(t) of the unique discrep¬ 
ancy between two coupled systems that initially differ by exactly one particle (see 
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[lq . part III). This makes sense for example for Hammersley’s process and exclu¬ 
sion type processes. If the particle system lives in one-dimensional space, we can 
also look at the system as the height function of an interface, so that the occupa¬ 
tion variables of the particle system define the increments of the height function. 
In terms of the height functions, we start by coupling two systems a and £ so that 
initially ( = a to the left of A'(0), and £ = a + 1 to the right of A(0). Then at all 
times the point X(t) is the boundary of the set {a; : ((x,t) = a(x,t) + 1}. 

This last idea generalizes naturally to the multidimensional interface model. We 
couple two height processes a and ( that satisfy a < ( < a -f 1 at all times, and 
prove that the boundary of the random set {a; : C(ie, t) = a(x,t) + 1} follows the 
characteristics of the macroscopic equation. 

Laws of large numbers for height functions of asymmetric interface models of 
the general type considered here have been earlier studied in a handful of papers. 
A hydrodynamic limit for ballistic deposition was proved in |24J, and for models 
that generalize the exclusion process in 3, 1 . 20 ]. Articles H, 24| deal with totally 
asymmetric models, and utilize the path-level variational formulation that general¬ 
izes from one-dimensional systems [23? . Article [2fJ introduces a different approach 
for partially asymmetric systems. These results are existence results only. In other 
words convergence to a limiting evolution is shown but nothing explicit about the 
limit is known, except that it is defined by a Hamilton-Jacobi equation. For partially 
asymmetric systems in more than one dimension it is presently not even known if 
the limit is deterministic. 

Our motivation for introducing a new model is to have a system for which better 
results could be proved. An advantage over earlier results is that here we can write 
down explicitly the macroscopic velocity function up to a constant factor. This 
is because the process is constructed through a homogeneous space-time Poisson 
process, so we can simultaneously scale space and time. This is not possible for 
a lattice model. With an (almost) explicit velocity we can calculate macroscopic 
profiles, for example see what profiles with shocks and rarefaction fans look like. In 
the earlier cases at best we know that the velocity function is convex (or concave), 
but whether the velocity is C 1 or strictly convex is a hard open question. Here this 
question is resolved immediately. 

Hammersley’s process has been a fruitful model for studying large scale behav¬ 
ior of one-dimensional asymmetric systems, by virtue of its connection with the 
increasing sequence model. For example, by a combination of the path-level varia¬ 
tional construction and the Baik-Deift-Johansson estimates jH, one can presently 


prove the sharpest out-of-equilibrium fluctuation results for this process [25), 127]. 


The model introduced in the present paper has a similar connection with a simple 
combinatorial model, so it may not be too unrealistic to expect some benefit from 
this in the future. 

Before the arrival of the powerful combinatorial and analytic approach pioneered 
in [3J, Hammersley’s process was used as a tool for investigating the increasing se¬ 
quences model. This approach was successful in finding the value of the limiting 
constant UM and in large deviations [22[ ■ The proofs relied on explicit knowledge 
of invariant distributions of Hammersley’s process. A similar motivation is possible 
for us too, and this time the object of interest would be the height of the ran¬ 
dom partial order. But currently we have no explicit knowledge of steady states of 
the process introduced here, so we cannot use the process to identify the limiting 
constant for the random order model. 

Recently Cator and Groeneboom || developed an approach to the one-dimensio¬ 
nal Hammersley process that captures the correct order t 1 / 3 of the current fluctua- 
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tions across a characteristic. The argument utilizes precise equilibrium calculations 
and a time reversal that connects maximizing increasing paths with trajectories 
of second class particles. In Q this method is adapted to the totally asymmetric 
exclusion process. Whether the idea can be applied in multidimensional settings 
remains to be seen. 

Organization of the paper. We begin by reminding the reader of the random 
partial order model, and then proceed to define the process and state the limit 
theorems. Proofs follow. A technical appendix at the end of the paper shows that 
the process has a natural state space that is a complete, separable metric space. 

2. The height of a random partial order 

Fix an integer v > 2. Coordinatewise partial orders on R" are defined for points 
x = (x!,...,x v ) and y = {y x ,..., y v ) by 

(1) x < y iff Xi < yi for 1 < j < v ) and x < y iff Xi < yi for 1 < j < u. 

We use interval notation to denote rectangles: (a, 6] = G R^ : a < x < b} for 
a < b in R 8 ', and similarly for [a, b] and the other types of intervals. 

Consider a homogeneous rate 1 Poisson point process in RA A sequence of Pois¬ 
son points p k , 1 < k < m, is increasing if p 1 < p 2 < ■ ■ ■ < p m in the coordinatewise 
sense. For a < b in R y , let H(a, b ) denote the maximal number of Poisson points on 
an increasing sequence contained in the set (a, b]. Let 1 = (1,1,..., 1) £ RA (This 
is the only vector we will denote by a boldface.) Kingman’s subadditive ergodic 
theorem and simple moment bounds imply the existence of constants c v such that 

(2) lim — H(0,nl) = c„ a.s. 

n—>oo yi 

Presently the only known value is = 2, first proved by Vershik and Kerov [29| 
and Logan and Shepp E3 . The case v = 2 is the same as the problem of the longest 
increasing subsequence of a random permutation, see Q for a review. Bollobas and 
Winkler |7| proved that c v —> e as v —> oo. 

The general case is called the ^-dimensional random partial order, and H(0,nl) 
is the height of the random partial order. The study of random partial orders was 
initiated by Winkler [30| . Here is an alternative construction of the random partial 
order on a fixed (rather than Poisson) number of elements. From the k\ linear orders 
on a set of k elements, choose v orders -<i, -< 2 , ■ ■ ■ , uniformly at random with 
replacement. Define the random order -< as the intersection, namely x -< y iff x ~<j y 
for j = 1,..., v. The height of the random order is the maximal size of a linearly 
ordered subset. Conditioned on the number k of Poisson points in the cube (0,nl], 
H(0, nl) has the same distribution as the height of the random order 

Let us also point out that by the spatial scaling of the Poisson point process, for 
any b = (bi ,..., b v ) > 0 in R d , 

(3) lim — H(0, nb) = Cj,(&i& 2&3 • • • b^) 1 ^ a.s. 

n —>00 YI 

3. The interface process 

Fix a spatial dimension d > 2. Appropriately interpreted, everything we say is 
true in d = 1 also, but does not offer anything significantly new. We describe 
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the evolution of a random, integer-valued height function a = (<7(2:)) lER d. Height 
values ±00 are permitted, so the range of the height function is Z* = Z U {± 00 }. 
The state space of the process is the space £ of functions a : R d —► Z* that satisfy 
conditions (i)-(iii): 

(4) (i) Monotonicity: x < y in R d implies cr(x) < cr(y). 

The partial order x < y on R d is the coordinatewise one defined in Section [2] 

(ii) Discontinuities restricted to a locally finite, countable collection of coordinate 
hyperplanes: for each bounded cube [—gl,gl] C R d , there are finite partitions 

—q = < si < ■ ■ ■ < s™* = q 

along each coordinate direction (1 < i < d), such that any discontinuity point of er 
in [— ql, gl] lies on one of the hyperplanes {x G [—gl, gl] : x, = s^}, 1 < i < d and 
0 < k < rrii. 

At discontinuities a is continuous from above: a(y) —> a(x) as y —> x so that 
y > x in R rf . Since a is Z*-valued, this is the same as saying that a is constant on 
the left closed, right open rectangles 


(5) 


k „/c+l 


11 -■ 


M+l 


), 


d 

k = (ki,k 2 ,...,kd) G JJ{0,1,2 ,...,rm - 1}, 
2—1 


determined by the partitions {.sf : 0 < k < rrii }, 1 < i < d. 

(iii) A decay condition “at — 00 ”: 

(6) for every b G R d , lim sup { \y\^^ d+1 ' 1 (x(y) : y < b, |y|oo > m\ = — 00 . 

M—> 00 l J 

The role of the (arbitrary) point b in condition ([6]) is to confine y so that as the 
limit is taken, all coordinates of y remain bounded above and at least one of them 
diverges to — 00 . Hence we can think of this as “y —> — 00 ” in R d . The £°° norm on 

is |y|oo = maxi<j< d \yi\. 

We can give £ a complete, separable metric. Start with a natural Skorohod 
metric suggested by condition (ii). On bounded rectangles, this has been considered 
earlier by Bickel and Wichura |5j, among others. This metric is then augmented with 
sufficient control of the left tail so that convergence in this metric preserves (0 • The 
Borel cr-field under this metric is generated by the coordinate projections a 1 —► <j{x). 
These matters are discussed in a technical appendix at the end of the paper. 

Assume given an initial height function a G £. To construct the dynamics, 
assume also given a space-time Poisson point process on R d x (0,oo). We define 
the process a(t) = {cr(x, t) : x G R d } for times t G [0, 00 ) by 

(7) a(x, t) = sup {< 7 ( 2 /) + H((y, 0), (x, t))}. 

y:y<x 

The random variable H((j/,0), (x,t)) is the maximal number of Poisson points on 
an increasing sequence in the space-time rectangle 

((y, 0), (x, f)] = {(77, s) G R d x (0, t]:yi< <Xi(l<i< d)}, 

as defined in Section [2j One can prove that, for almost every realization of the 
Poisson point process, the supremum in © is achieved at some y, and a(t) G £ for 
all t > 0. In particular, if initially cr{x) is finite then a(x,t) remains finite for all 
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0 < t < oo. And if a(x) = ±oo, then a(x,t) = a(x) for all 0 < t < oo. This defines 
a Markov process on the path space Z?([0, oo), E). 

The local effect of the dynamical rule 0 is the following. Suppose (y, t ) £ 
R d x (0, oo) is a Poisson point, and the state at time t— is er. Then at time t the 
state changes to cr y defined by 

/ox _!//>, _ / <r(x) + l,Hx>y and a(x) = a(y), 

' ' ' ]_cr(x), for all other x £ R d . 


We can express the dynamics succinctly like this: Independently at all x £ R d , 
cr(x) jumps to cr(x) + 1 at rate dx (d-dimensional volume element). When a jump 
at x happens, the height function a is updated to er + 1 on the set {u> £ R' ( : w > 
x,cr(w ) = cr{x)} to preserve the monotonicity property 0. It also follows that if 
cr(y) = ±oo then o v = a. 

We express this in generator language as follows. Suppose 0 is a bounded mea¬ 
surable function on E, and supported on a compact cube K C R d . By this we mean 
that 0 is a measurable function of the coordinates (ct(x)) x ek- Define the generator 
£ by 

(9) Ccj)(a) = f [<j>(a y ) - <j>(cr)]dy. 

J R d 

The next theorem verifies that C gives the infinitesimal description of the process 
in one basic sense. 

Theorem 3.1. For bounded measurable functions <f> on E, er £ E, and t > 0, 

(10) E a [(f>(<j(t))} — 4>(a) = f E a [C(/)(a{s))}ds. 

Jo 

E a denotes expectation under the path measure P a of the process defined by 0 and 
started from state a. 


4. Hydrodynamic limit for the height process 

Let it 0 : R d —> R be a nondecreasing locally Lipschitz continuous function, such 
that for any b £ R d , 

(11) lira sup{|y|^ d/(d+1) u 0 ( 2 /) : y < b, \y\ x > m\ = -oo. 

M —»oo k ) 

The function uq represents the initial macroscopic height function. Assume that 
on some probability space we have a sequence of random initial height functions 
{cr n (y, 0) : y £ R d }, indexed by n. Each a n (■ , 0) is a.s. an element of the state space 
E. The sequence satisfies a law of large numbers: 

(12) for every y £ R d , n~ 1 a n (ny, 0) —> uo(y) as n —* oo, a.s. 

Additionally there is the following uniform bound on the decay at — oo: 

for every fixed b £ R d and C > 0, with probability 1 there exist finite, 

(13) possibly random, M,N > 0 such that, if n > N, y < b and |y|oo > M, 
then a n {ny, 0) < -Cn\y\ d JJ' d+1) ■ 
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Augment the probability space of the initial cr ra (-,0) by a space-time Poisson 
point process, and define the processes a n (x,t) by 0. For x = (aq,... ,Xd) > 0 in 
R d , define 

g(x) = c d +i{xix 2 x 3 • • ■x d ) 1,lyd+l) ■ 

The constant Cd+i is the one from (|2|) . and it comes from the partial order among 
Poisson points in d + 1 dimensional space-time rectangles. 

Define a function u(x,t) on R d x [0, oo) by u(x,0) = u o(x) and for t > 0, 

(14) u(x, t ) = sup {u 0 (y) + tg{( x - y)/t)}. 

y.y<x 

The function u is nondecreasing in x, increasing in t , and locally Lipschitz in R d x 
(0, oo). 

Theorem 4.1. Suppose uq is a locally Lipschitz function on R d that satisfies CD- 
Define u(x,t) through (11411 . Assume that the initial random interfaces Wn{y, 0)} 
satisfy (fT21) and (fldl) . Then for all ( x,t ) £ R d x [0, oo), 

(15) lim n~ 1 a n (nx,nt) = u(x,t) a.s. 


By the monotonicity of the random height and the continuity of the limiting 
function, the limit (fl5l) holds simultaneously for all (x, t) outside a single exceptional 
null event. 

Extend g to an u.s.c. concave function on all of R d by setting g = —oo outside 
[0, oo) d . Define the constant 


(16) 


Hd 


Cd+l 

d + 1 


d+1 


The concave conjugate of g is g*{p) = inf^ja; • p — g{x)}, p £ R d . Let / = — g*. 
Then f{p) = oo for p £ (0, oo) d , and 

(17) f(p) = n d (p 1 p 2 ■ ■ ■ Pd)~ l for p > 0 in R d . 


The Hopf-Lax formula (flTll implies that u solves the Hamilton-Jacobi equation (see 

0 ) 


(18) 


d t u - /(Vm) = 0, u| t= o=it 0 . 


In other words, /(Vu) is the upward velocity of the interface, determined by the 
local slope. 

The most basic case of the hydrodynamic limit starts with er(y, 0) = 0 for y > 0 
and a(y, 0) = — oo otherwise. Then a(x,t) = H((0,0), (x, t)) for x > 0 and —oo 
otherwise. The limit is u(x,t) = tg{x/t). 


5. The defect boundary limit 

Our objective is to generalize the notion of a second class particle from the one¬ 
dimensional context. The particle interpretation does not make sense now. But a 
second class particle also represents a defect in an interface, and is sometimes called 
a ‘defect tracer.’ This point of view we adopt. Given an initial height function 
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<x(y, 0), perturb it by increasing the height to a(y, 0) + 1 for points y in some set 
A(0). The boundary of the set A(0) corresponds to a second class particle, so we 
call it the defect boundary. How does the perturbation set A(-) evolve in time? To 
describe the behavior of this set under hydrodynamic scaling, we need to look at 
how the Hamilton-Jacobi equation (1181) carries information in time. 

For (x,t) £ R d x (0, oo), let I(x,t) be the set of maximizers in (fill) : 

(19) = {y £ R d : y < x, u(x,t) = u 0 {y) +tg({x - y)/t)}. 

Continuity and hypothesis m guarantee that I(x,t) is a nonempty compact set. 
It turns out that these three statements (i)-(iii) are equivalent for a point ( x,t ): 
(i) the gradient Vw in the a;-variable exists at (x,t), (ii) u is differentiable at (x,t), 
and (iii) I(x,t) is a singleton. We call a point (x,t) with t > 0 a shock ii I(x,t) has 
more than one point. 

For y £ R d let W (y, t ) be the set of points x £ R d for which y is a maximizer in 
the Hopf-Lax formula (fT4l) at time t: 

(20) W(y,t) = {x £ R d : x > y, u(x,t ) = u 0 (y) + tg((x - y)/t)}, 
and for any subset B C R d , 

(21) W(B,t) = U W(y,t). 

y£B 


Given a closed set B C R rf , let 


( 22 ) X(B,t) = W{B,t) nW(B c ,t). 

W(B,t ) and W(B c ,t) are both closed sets. We can characterize x £ X(B,t) as 
follows: if (x,t) is not a shock then the unique maximizer {y} = I(x,t) in (|14[l lies 
on the boundary of B, while if (x,t) is a shock then I(x,t) intersects both B and 

w. 

If dimension d = 1 and B = [a, oo) C R, an infinite interval, then X(B,t) is 
precisely the set of points x for which there exists a forward characteristic x(-) 
such that a;(0) = a and x(t) = x. By a forward characteristic we mean a Filippov 
solution of dx/dt = f'(S7u(x , t)) [9(,ll8|. A corresponding characterization ofX(B, t ) 
in multiple dimensions does not seem to exist at the moment. 

The open e-neighborhood of a set B C R d is denoted by 

(23) B^ = {x : d(x, y) < e for some y £ B}. 


The distance d(x, y) can be the standard Euclidean distance or another equivalent 
metric, it makes no difference. Let us write B^~ e ^ for the set of x £ B that are at 
least distance e > 0 away from the boundary: 


(24) 


^ = {x £ B : d(x, y) > £ for all y £ B} 


(H c )( £ ) 


The topological boundary of a closed set B is bdR = B fl B c . 

Suppose two height processes a(t ) and £(i) are coupled through the space-time 
Poisson point process. This means that on some probability space are defined the 
initial height functions cr(y, 0) and £(?/, 0), and a space-time Poisson point process 
which defines all the random variables H((y, 0), (x, t)). Process cr(x,t) is defined 
by 0, and process ((x,t) by the same formula with a replaced by £, but with 
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the same realization of the variables H((y, 0), ( x , t)). If initially a < £ < a + h for 
some constant h, then the evolution preserves these inequalities. We can follow the 
evolution of the “defect set” A(t), defined as A(t) = {x : £( x,t ) = a(x,t ) + h} for 
t > 0. This type of a setting we now study in the hydrodynamic context. In the 
introduction we only discussed the case h = 1, but the proof works for general finite 
h. 

Now precise assumptions. On some probability space are defined two sequences of 
initial height functions a n (y,0) and C, n (y, 0). The {a n {y, 0)} satisfy the hypotheses 
m and [THill of Theorem 14.11 For some fixed positive integer h , 

(25) <r n (y,0) < ( n (y, 0) < a n (y, 0) + h for all n and y G R rf . 

Construct the processes a n (t) and {t) with the same realizations of the space-time 
Poisson point process. Then 

(26) cr n (x,t) < ( n (x,t) < cr n (x,t) + h for all n and (x,t). 

In particular, and a n satisfy the same hydrodynamic limit. 

Let 

(27) A n {t) = {ieR i : ( n (x,t) = a n (x,t ) + h}. 

Our objective is to follow the evolution of the set A n (t) and its boundary bd{A„(f)}. 
We need an initial assumption at time t = 0. Fix a deterministic closed set B C R d . 
We assume that for large n, n~ 1 A n { 0) approximates B locally, in the following sense: 
almost surely, for every compact K C R d and e > 0, 

(28) B O K C {n -1 A n (0)} OA' C B^ fl K for all large enough n. 


Theorem 5.1. Let again uq satisfy (ED and the processes a n satisfy m and (USD 
at time zero. Fix a positive integer h and a closed set B C R t? . Assume that the 
processes o n are coupled with processes through a common space-time Poisson 
point process so that (12611 holds. Define A n (t) by (12711 and assume Al n (0) satisfies 

U2SD. 

IfW(B,t) = 0, then almost surely, for every compact I\ C R d , A n (nt)DnK = 0 
for all large enough n. 

Suppose W(B,t) if 0- Then almost surely, for every compact K C R d and e > 0, 

(29) bd {?r _1 A n (nt)} fl K C X(B, t)^ 0 K for all large enough n. 

In addition, suppose no point ofW(B c ,t) is an interior point ofW(B,t). Then 
almost surely, for every compact K C R d and e > 0, 

W{B,t)^ n K C {n^A^nt)} n K 

C W(B,t)^ O K for all large enough n. 


The additional hypothesis for (1301) . that no point of W(B C , t) is an interior point 
of W(B,t), prevents B and B c from becoming too entangled at later times. For 
example, it prohibits the existence of a point y £ bdB such that W(y,t ) has 
nonempty interior (“a rarefaction fan with interior”). 
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6. Examples and technical comments 


6.1. Second class particle analogy 


Consider a one-dimensional Hammersley process z(t) = (z»(£))iez with labeled 
particle locations • ■ ■ < Z-±(t) < zo(t) < zi(£) < ■■■ on R. In terms of labeled 
particles, the infinitesimal jump rule is this: Zj jumps to the left at exponential rate 
Zi — Zi—i, and when it jumps, its new position z\ is chosen uniformly at random from 
the interval (zi- 1 , Zi). The height function is defined by <r(x, t ) = sup{i : zft) < x} 
for itR. 

Now consider another Hammersley process 5(f) constructed with the same re¬ 
alization of the space-time Poisson point process as z(t). Assume that at time 0, 
5(0) has exactly the same particle locations as z(0), plus h additional particles. 
Then at all later times 5(f) will have h particles more than z(t), and relative to the 
z(£)-process, these extra particles behave like second class particles. 

Suppose the labeling of the particles is such that 5j(£) = zft) to the left of all 
the second class particles. Let X- t (£) < • • • < Xh(t) be the locations of the second 
class particles. Then the height functions satisfy a(x,t) = er(x,£) for x < X\ (f), and 
ct(x, t ) = a(x, t ) + h for x > Xh{t). So in this one-dimensional second class particle 
picture, the set A(t) is the interval [Xh(t), oo). It has been proved, in the context of 
one-dimensional asymmetric exclusion, A’-exclusion and zero-range processes, that 
in the hydrodynamic limit a second-class particle converges to a characteristic or 
shock of the macroscopic p.d.e. 


12, 18 


Despite this analogy, good properties of the one-dimensional situation are readily 
lost as we move to higher dimensions. For example, we can begin with a set A(0) 
that is monotone in the sense that x £ A(0) implies y £ A(0) for all y > x. 
But this property can be immediately lost: Suppose a jump happens at w such 
that C( w >0) = a(w,0) but the set V = {x > w : a(x, 0) = cr(w, 0)} intersects 
A(0) = {x : C( x >0) = ct(x, 0) + 1}. Then after this event C, = a on V, and cutting 
V away from A(0) may have broken its monotonicity. 


6.2. Examples of the limit in Theorem 1 5. 1\ 


We consider here the simplest macroscopic profiles for which we can explicitly calcu¬ 
late the evolution W(B , t) of a set B , and thereby we know the limit of n -1 A n (nt) in 
Theorem o These are the flat profile with constant slope, and the cases of shocks 
and rarefaction fans that have two different slopes. Recall the slope-dependent 
velocity f(p) = Kd(piP 2 P 3 • • ■ Pd)^ 1 for p £ (0, oo) d , where Kd is the (unknown) 
constant defined by (J2]) and (USD- 

For the second class particle in one-dimensional asymmetric exclusion, these 
cases were studied in mm. 

Flat profile. Fix a vector p £ (0, oo) d , and consider the initial profile uq(x) = 
p ■ x. Then u(x,t) = p ■ x + tf(p), for each (x, f) there is a unique maximizer 
y(x,t) = x + tVf(p) in the Hopf-Lax formula, and consequently for any set B , 
W(B,t) = -tXf(p) + B. 

Shock profile. Fix two vectors A, p £ (0, oo) d , and let 


(31) 


ito(x) 


J p ■ x, (p — A) • x > 0, 
( A ■ x, (p — A) ■ x < 0. 
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Then at later times we have 
u(x,t) = 

The Hopf-Lax formula is maximized by 


p-x + tf(p), (p - X)-x> t(f( A) - f(p)), 
X-x + tf( A), (p - A) • x < t(f{ A) - /(p)). 


jx + tVf(p), if (p - A) • x > t(f(X) - f(p)), 
1 x + tV/( A), if (p-X)-x< t(/( A) - /(p)). 


In particular, points (x, t) on the hyperplane (p — A) • x = t(f(X) — f(p)) are shocks, 
and for them both alternatives above are maximizers. In the forward evolution, 
W ( y , t) is either a singleton or empty: 


' y-tvf(p), 

if (p — A) • y > t(f( A) - /(p)) + i(p - A) • V/(p), 

0, if t(f(X) - f{p)) + t(p - A) • V/(A) < (p — A) ■ y 

< *(/(A) — /(p)) + t(p - A) • V/(p), 
A), 

if (P - A) • y < - /(p)) + t(p - A) • V/(A). 


In this situation Theorem 15.II is valid for all sets B. 

Rarefaction fan profile. Fix two vectors A, p £ (0, oo) d , and let 


uoO) 


J A • x, (p — A) • x > 0, 
[ p ■ x, (p — X) • x < 0. 


For (x, t) such that 

-t(p - A) • V/(p) < (p - A) • x < -t(p - A) • V/(A) 
there exists a unique s = s(x, f) £ (0,1) such that 

(p- X) ■ x = -t(p - A) • V/(sA + (1 - s)p). 


Then at later times the profile can be expressed as 

f P-X + tf(p), if (p - A) • X < -t(p - A) • V/(p), 

,./ r ^ = J (sA + (1 - s)p) • a; + t/(sA + (1 - s)p), if 
l , j 1 ~t(p ~ A) • V/(p) < (p — X) ■ x < —t(p — A) • V/(A), 

[ A • x + tf( A), if (p - X)-x> -t(p - A) • V/(A). 

The forward evolution manifests the rarefaction fan: points p on the hyperplane 
(p — A) • y = 0 have W ( y , t) given by a curve, while for other points y W(y. t ) is a 
singleton: 

f y-*V/(p), if (p- A) -y < 0, 

^(j/> ^ = ] {p — tV/(sA + (1 — s)p) : 0 < s < 1}, if (p — A) • y = 0, 
l y - iV/(A), if (p - A) • y > 0. 

In Theorem 15.II consider the half-space B = {x : (p — A) • x > 0}. Then 
X{B,t) = {x : -t{p - A) • V/(p) <(p — X) - x < -t{p - A) • V/(A)}, 
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the “rarefaction strip” in space. Statement (®D is not valid for B, because the 
interior of X(B,t) lies in the interiors of both W(B,t) and W(B c ,t). Statement 
(l29l) is valid, and says that the boundary of n~ 1 A n (nt) is locally contained in any 
neighborhood of X(B,t). 

In the corresponding one-dimensional setting, Ferrari and Kipnis j~K^ proved that 
on the macroscopic scale, the second class particle is uniformly distributed in the 
rarefaction fan. Their proof depended on explicit calculations with Bernoulli distri¬ 
butions, so presently we cannot approach such precise knowledge of bd{n” 1 2 l rl (nt)}. 


6.3. Some random initial conditions 

We give here some natural examples of random initial conditions for Theorems 14.11 
and lb.ll for the case d = 2. We construct these examples from space-time evolutions 
of one-dimensional Hammersley’s process. The space-time coordinates (y, t) of the 
1-dimensional process will equal the 2-dimensional spatial coordinates x = (x±,X 2 ) 
of a height function. 

Flat profiles. In one dimension, Aldous and Diaconis [ij denoted the Hammersley 
process by N(y,t). The function y e-> N(y,t) (y € R) can be regarded as the 
counting function of a point process on R. Homogeneous Poisson point processes 
are invariant for this process. 

To construct all flat initial profiles Uq(x) = p ■ x on R 2 , we need two parameters 
that can be adjusted. The rate p of the spatial equilibrium of N(y,t) gives one 
parameter. Another parameter r is the jump rate, in other words the rate of the 
space-time Poisson point process in the graphical construction of N(y,t). Let now 
N(y,t) be a process in equilibrium, defined for —oo <t< oo, normalized so that 
TV(0, 0) =0, with jump rate r, and so that the spatial distribution at each fixed time 
is a homogeneous Poisson process at rate p. Then the process of particles jumping 
past a fixed point in space is Poisson at rate r/p [l|, Lemma 8. Consequently 
EN{y,t) = py+(r/p)t. 

This way we can construct a random initial profile whose mean is a given flat 
initial profile: given p = (phpfi) £ (0,oo) 2 , take an equilibrium process {N(y,t) : 
y G R,f G R} with p = pi and r = P 1 P 2 , and define the initial height function for 
x = {xi,xfi) £ R 2 by cr((xi,X 2 ), 0) = N(xi,x 2 ). 

Shock profiles. Next we construct a class of initial shock profiles. Suppose p = 
{pi , P 2 ) and A = (Ai, A 2 ) satisfy p > A and pi/ P 2 < A 1 /A 2 . Start by constructing the 
equilibrium Hammersley system {N(y, t) : y E R, t £ R} with spatial density p = 
Ai and jump rate r = A 1 A 2 . Set a = (p\ — \i)/(p 2 — A 2 ) > 0. Stop each Hammersley 
particle the first time it hits the space-time line t = —ay, and “erase” the entire 
evolution of N(y, t) above this line. The assumption pi/P 2 < A 1 /A 2 guarantees that 
each particle eventually hits this line. Now we have constructed the slope-A height 
function cr((xi, X 2 ), 0) = N(x 1 , 2 : 2 ) below the line (p — A) ■ x = 0 X 2 = —ax 1 . 
(Slope-A in the sense that Ea(x, 0) = A • x.) 

To continue the construction, put a rate t' = P 1 P 2 space-time Poisson point 
process above the line t = —ay in the space-time picture of the 1-dim Hammersley 
process. Let the Hammersley particles evolve from their stopped locations on the 
line t = —ay, according to the usual graphical construction [l| of the process, using 
the rate r' space-time Poisson points. The construction is well defined, because 
given any finite time T, N(y,T) is already constructed for y < —T/a, and for 
y > —Tla the particle trajectories can be constructed one at a time from left to 
right, starting with the leftmost particle stopped at a point (y, —ay) for y > —T/a. 
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One can check that defining cr((x 1 ,^ 2 ),0) = N(x 1 , 2 : 2 ) for X 2 > —cix 1 gives the 
slope-p height function above the line (p — A) • x = 0. Now we have a random initial 
height function er(x,0) with mean Ea(x,0) = uq(x) as in (IdTl) . 

Finally, we describe a way to define initial configurations for the coupled pro¬ 
cesses ( and a in the context of this shock example. We shall do it so that the set 
{x : (C(x,0) = <r(x,0) + 1} lies inside B = {x : X 2 > —ax 1 }, and approximates 
it closely. Let 0) be the height function defined above in terms of the N(y,t) 
constructed in two steps, first below and then above the line t = —ay. Let Zk(t) 
be the trajectories of the labeled Hammersley particles. These trajectories are the 
level curves of £(x,0), namely £((xi, X 2 ), 0) > k iff Zk(x 2 ) < x\. The construction 
performed above has the property that each Zk(t) crosses the line t = —ay exactly 
once (the particles were stopped upon first hitting this line, and then continued 
entirely above the line). 

Define new trajectories z' k (t ) as follows: z k (t) = Zk(t) below the line t = —ay. 
From the line t = —ay the trajectory z k (t) proceeds vertically upward (in the t- 
direction) until it hits the trajectory of Zk+i(t). From that point onwards z k (t) 
follows the trajectory of z k +\ (t). This is done for all k. Let N'(y, t) be the counting 
function defined by N'(y,t) = supjfc : z k (t) < y}. And then set ct((xi, X2), 0) = 
N'(x 1 ,x 2 ) 

The initial height functions er(x, 0) and £( x , 0) thus defined have these properties: 
< 7 ( 2 ;, 0) = £(x,0) f° r %2 < —ax 1 . For any point (xi,X 2 ) such that 2:2 > — 02:1 and 
some particle trajectory Zk(t) passes between (xi,— 02 : 1 ) and (xi,X 2 ), C( a; ?0) = 
cr(x,0) + 1. This construction satisfies the hypotheses of Theorem 15.II 

6 . 4 ■ Some properties of the multidimensional Hamilton-Jacobi 
equation 

Let u(x,t) be the viscosity solution of the equation ut = f(Vu), defined by the 
Hopf-Lax formula (1141) . By assumption, the initial profile uq is locally Lipschitz 
and satisfies the decay estimate CD- Hypothesis m is tailored to this particular 
velocity function, and needs to be changed if / is changed. 

Part (b) of this lemma will be needed in the proof of Thm. 15.11 

Lemma 6.1. (a) For any compact K C R d ; \J x K I(x,t) is compact. 

(b) W(B,t) is closed for any closed set B C R d . 

Proof, (a) By (fill) , as y —> —00 for y < x, uo(y) + tg((x — y)/t) tends to —00 
uniformly over x in a bounded set. Also, the condition inside (1191) is preserved by 
limits because all the functions are continuous, (b) If W(B 1 t) 9 Xj —> x, then by 
(a) any sequence of maximizers yj € I (xj ,t)r\B has a convergent subsequence. □ 

The association of I(x,t) to x is not as well-behaved as in one dimension. For 
example, not only is there no monotonicity, but a simple example can have xi < X 2 
with maximizers yi € I ( 2 , t) such that 2/2 < yi ■ The local Lipschitz condition on 
uq guarantees that each y G I{x,t) satisfies y < x (i.e. strict inequality for all 
coordinates). 

Properties that are not hard to check include the following. Part (a) of the lemma 
implies that u(x, t ) is locally Lipschitz on R d x (0, 00 ). Lipschitz continuity does not 
necessarily hold down to t = 0, but continuity does, u is differentiable at (x, t) iff 
J(x,t) is a singleton {y}, and then Vtt(x,t) = V</((x — y)/t). Also, Vu is continuous 
on the set where it is defined because whenever (x n , t n ) —> (x, t ) and y n £ I(x n ,t n ), 
the sequence {y n } is bounded and all limit points lie in I(x,t). 
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A converse question is when W(y,t) has more than one point. As in one dimen¬ 
sion, one can give a criterion based on the regularity of u o at y. The subdifferential 
D~uq(x) and superdifferential D + uo(x ) of uq at x are defined by 

D-M*) = 1« 6 : lim tot "°fa> - “»(») - g ■ (y - *> > „l 

l v-x \\y-x\\ J 

and 

B+«„(x) = Le‘: lim sup "°(ri ~ - P ■ (v ~ *) £ „) 

l lly — X ll J 

It is a fact that both D±uq(x) are nonempty iff ito is differentiable at x , and then 
D ± u 0 (x) = {VmoOO}- 

One can check that W(y,t) C y — tV/(Z) + u 0 (j/)). Consequently if D~u 0 (y) 
is nonempty, W(y,t) cannot have more than 1 point. Another fact from one¬ 
dimensional systems that also holds in multiple dimensions is that if we restart 
the evolution at time s > 0, then all forward sets W(y,t) are empty or singletons. 
In other words, if u is a solution with initial profile u o, and we define uq(x ) = u(x, s) 
and u(x, t ) = u(x, s+t ), then D~uo(y) is never empty. This is because 'Vg((x — y)/s) 
lies in s)}(x) for every y that maximizes the Hopf-Lax formula for u(x, s). 


7. Proof of the generator relation 


In this section we prove Theorem 13.II Throughout the proofs we use the abbrevia¬ 
tion 

X\ = X1X2X3 ■ ■ ■ Xd 

for a point x = (x\, ..., Xd) G R d . We make the following definition related to the 
dynamics of the process. For a height function a e £ and a point x 6 R d , let 


(32) 


S x (<r) 


{y £ R d : y < x,cr(y) = cr(x)}, if cr(x) is finite, 
0 , if <j(x) = ± 00 . 


S x (a) is the set in space where a Poisson point must arrive in the next instant 
in order to increase the height value at x. Consequently the Lebesgue measure 
(volume) | S x (cr) | is the instantaneous rate at which the height a(x ) jumps up by 
1. Since values a(x) = ±00 are not changed by the dynamics, it is sensible to set 
S x (a) empty in this case. For a set K in R d we define 


(33) 


Sk{<t) = (J S x (a), 

xG K 


the set in space where an instantaneous Poisson arrival would change the function 
a in the set K. 

We begin with a simple estimate. 

Lemma 7.1. Let x > 0 in R d , t > 0, and k a positive integer. Then 
P{H((0,0), (x, t)) > *} < dA)A < 

where the second inequality is valid if k > e 2 (xH) 1 ^ d+1 \ Note that above (0,0) 
means the space-time point (0, 0) € R“ x [0, 00 ). 
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Proof. Let 7 = x\t = X 1 X 2 X 3 ■ ■ ■ Xd • t > 0 be the volume of the space-time rectangle 
(0, a;] x (0, t]. k uniform points in this (d+l)-dimensional rectangle form an increasing 
chain with probability ( k\)~ d . Thus 

-P{H((0,0), (x, t)) > k} < Y, e -^f( J \k\)- d = 1 k {k\)-^ 

j:j>k J ' ' ' 

< 7 k (k/e)~ Hd+1) < e ~ k{d+1) 

if k > e 2 7 1/(d+1) . □ 

We need to make a number of definitions that enable us to control the height 
functions treE. For b £ R d and h £ Z, let y b ’ h (a) be the maximal point y < b such 
that the rectangle [y, 6] contains the set {x < b : cr( x) > h}, with y b,h (a ) = b if this 
set is empty. Note that if y b,h (a ) ^ b then there exists x < b such that a(x) > h 
and | b - x ^ = \b - y b ' h (a)\ OQ . 

Throughout this section we consider compact cubes of the type 

K = [~ql,ql] C R rf 

for a fixed number q > 0. When the context is clear we may abbreviate y h = 
y ql ’ h (cr). Define 

A fc (cr) = sup (ql-y h )\-(k-h)-( d+1 l 

— oo<h<.k —2 

Property (JG]) of the state space guarantees that Afe(cr) < 00 . 

The minimal and maximal finite height values in K are defined by 

I{K , cr) = min{cr(a;) : x £ K, —00 < cr(x) < 00 } 

and 

J(K,a) = max{cr(x) : x £ K, —00 < a(x) < 00 }. 

If cr = ±00 on all of K we interpret I(K,a) = 00 = —J(K,cr). Otherwise these 
quantities are finite because a can take only finitely many values in K. If a is finite 
on all of K then by monotonicity I(K,cr) = cr(—ql) and J(K,a) = a(ql). Set 

( J(K,a)+l 
fe=/(if,cr)+l 

If cr = ±00 on all of K then ^k(o') = 0. 

The next two lemmas are preliminary and illustrate how iI>k(<x) appears as a 
bound. 



Lemma 7.2. For a cube K = [— ql, ql] and a £ £, 
(35) |Sif(o-)| < 2 d+1 i> K (a). 


Proof. If cr = ±00 on all of K then both sides of (l35l) are zero by the definitions. 
Suppose I{K , cr) is finite (this is the complementary case). If a : £ K and y £ S x (a ), 
then y < x < ql and cr(y) = cr(x) > I(K,a), and consequently y £ [y I( - K,a \ql]. 
This is true for an arbitrary point y £ Sk{<x). Since < y / ( K . cr ) j W e can 

weaken the conclusion to S K {o) C [y 1 ^- 1 ^ 1] to get 


|^(a)|<(gl-y™- 1 )! 


2 d+1 (ql -yl(K.a)- 1)1 

{I{K,a) + l-{I(K,a)-l)) d +^ 


< 2 d+1 \ I{Kt<T)+1 (a) < 2 d+1 Tjj K (a). 


□ 
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Lemma 7.3. Define the event 

G = {there exist x £ K and y £ R‘ ; such that y < x, 
(36) — oo < a(y) < a(x) — 1 < oo, 

and H ((y, 0), (a;, t)) > a(x) + 1 - a(y)}. 

Then for 0 < t < 1/(2 e d+1 ipK{<x)), 

P°(G) < 2 e 2{d+1) 4 2 K (a)t 2 . 


Proof. Let the index k run through the finite values of cr(x) + l in K , and h represent 
a(y). Then 

J{K,c r)+l 

P a {G) < Y. E P{my h ,0),(ql,t))>k-h}. 

k=I(K,c r) + l h<k -2 


By Lemma 17. II and the inequality j! > (j/e) 3 , for a fixed k the inner sum becomes 


E 

h<k -2 


((gl -/)!*) 

((fc - /i )!) d+1 


k—h 


^E 

J>2 




(iO 


d+l 


< 2(e d+i Afe(cr)t)^. 


The assumption on t was used to sum the geometric series. Now sum over k. □ 
We get the first bound on the evolution. 

Lemma 7.4. Let <j> be a bounded measurable function on E, supported on a com¬ 
pact cube K = [—gl,gl] C K d . This means that 4>(cr) depends on a only through 
(&{x))xeK■ Then there is a finite constant C = C( H^Hoo) such that, for all a £ E 
and t > 0, the quantity 


A t (a) = E a [4>(<j(t))] - 4 >{<t) - tC(f>{<j) 

satisfies the bound | A* (cr) | < Ct 2 ifff i {a), provided 0 < t < l/(2e d+1 xfK (cr)). 

Proof. We may assume that a is not ±oo on all of K. For otherwise from the 
definitions 

E a [<f>(a(t))\ = 4>{a) and C<t>(a) = 0, 

and the lemma is trivially satisfied. 

Observe that on the complement G c of the event defined in (1551) . 

a(x,t)= sup (cr(y) + H((y,0),(x,t))} 

J/6S x (o-)U{a:} 

for all x £ K. (The singleton {cc} is added to S x (a) only to accommodate those 
x for which cr(x) = ±oo and S x (a) was defined to be empty.) Consequently on 
the event G c the value cf(a(t)) is determined by a and the Poisson points in the 
space-time region Sk(&) x (0, £]. Let D :] be the event that Sk(ct) x (0 ,t\ contains 
j space-time Poisson points, and D 2 = (Dq U D\) c the event that this set contains 
at least 2 Poisson points. On the event Di , let Y £ R d denote the space coordinate 
of the unique Poisson point, uniformly distributed on ^(c). Then 

E a [Ha(t))] = 4(a) ■ P a (G c n D 0 ) + E[4(a Y ) • l G e nI3l ] + O ( P a (G ) + P a (D 2 )) 

= 4(a) + tC4(a) +t 2 O (4 k(v) + |Sic(c)| 2 ) ■ 
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To get the second equality above, use 

P a {Dj) = {j\)~ 1 {t\S K {a)\y exp(-t|5' i< -(cr)|), 

Lemma 17751 for bounding P a {G ), and hide the constant from Lemma 17751 and ||</>||oo 
in the O-terms. Proof of the lemma is completed by □ 

We insert here an intermediate bound on the height H. It is a consequence of 
Lemma 17. II and a discretization of space. 

Lemma 7.5. Fix t > 0, a £ (0,1/2), and f3 > e 2 i 1 /( d + 1 ). Then there are finite 
positive constants do, C\ and Ci such that, for any 9 > 9 q, 

P{ there exist y < x such that \y\oo >9, la; — y |oo > ce\y\oo, 

(37) 

and H((y,0), (x,t)) > (3\x - y |^ (d+1) } < C\ exp(-C' 2 6» d/(<i+1) ). 


For positive to, define 

(38) a m (x,t)= sup {a{y) + H{{y,0),{x,t))} 

y<x 

\y\oa<m 

Corollary 7.6. Fix a compact cube K C R rf , 0 < T < oo, and initial state 
a £ E. Then there exists a finite random variable M such that, almost surely, 
a{x,t) = cr M {x,t) for (x,f) £ K x [0,T]. 

Proof. If a(x ) = ±oo then y = a; is the only maximizer needed in the variational 
formula 0. Thus we may assume that I{K , cr) is finite. 

Fix a € (0,1/2) and (3 > e 2 T 1 /( d+1 ). By the boundedness of K, (1371) . and Borel- 
Cantelli there is a finite random M such that 


U((y,0) 1 (x,T))<^\x-y\ d J 1 ^ 

whenever x £ K and |a/1^ > M. Increase M further so that M > 1 + |ar| for all 
x £ K , and a(y) < —{2(3 + \I(K,cr)\ + l)\y\((^^ d+1 ' > for all y such that y < a; for 
some x £ K and \y\ao > M. 

Now suppose y < x, x £ K , a(x) is finite, and Ij/Ioq > M. Then 

a{y) + H {{y, 0), (x, t)) < -{2(3 + \I{K, a)\ + l)\y\^ d/{d+1) + (3\x - y \ d ^ d+1 ^ 

< — \I{K,cr)\ — 1 < cr(x) — 1 < a{x,t) — 1. 


We see that y cannot participate in the supremum in 0 for any (x,t) £ I\ x 

[0, T\. □ 

To derive the generator formula we need to control the error in Lemma 17.41 
uniformly over time, in the form A r (cr(s)) with 0 < s < t and a small r > 0. For 
a fixed k, Afc(cr(s)) is nondecreasing in s because each coordinate of y h decreases 
over time. For q > 0 and k £ Z introduce the function 


(39) 


^q,k{o) 


sup - 

x<ql (l V {k 


d 

oo 




d +1 ' 


A calculation that begins with 


Afc(cr) < g d /2+ sup 
h<.k —2 
y ql ’ h ^qi 


2 d \y ql ’ h \L 

{k - h) d + l 
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shows that 

A fe (a) <q d /2 + 2 d ^ qik (a). 

Interface heights cr(x, s ) never decrease with time, and ^ q ^{cr) is nonincreasing in k 
but nondecreasing in a. Therefore we can bound as follows, uniformly over s G [0, t]: 

J(K,c t(s)+1 

A fc( a ( s )) 

k=I (.ftr,f7(.s))-|-l 

< (j(K,cr(t)) — I(K,cr( 0)) + l) • max Xi(a(s)) 

(40) < (J(K, cr(t)) — I(K, cr(0)) + l) 2 (g 2d + 1) + 2 4d 4'^ J(if o . (0)) (cr(i)). 

Above we used the inequality c(a + b) 2 < 2ca 2 + c 2 + b 4 for a,b,c > 0. The next 
lemma implies that the moments E a [' i p q j, R .(cr(t))] are finite for all p < 00 . 

Lemma 7.7. Let a be an element of the state space E. Fix t > 0 and a point 
ql G R+. Then there exists a finite number vo(cr) such that, for v > vq( cr), 

(41) P a {'fr q j( K:a) (a(t)) >v}< C 1 exp(-C 2 v 1/( - d+1) ), 

where the finite positive constants C \, C 2 are the same as in Lemma | 7. r >\ above. 
Proof. Choose a , /3 so that (1371) is valid. Let 

(3\ = 2f3 + f3{2a) d/(d+1) + 2. 

Fix vq = vo(cr) > 0 so that these requirements are met: vq > 1 + | I(K, er)| d+1 , and 
for all y < x < ql such that |x|(^. > vq, 

a(x) < -(3 1 \x\ d 0 ^ d+1) and \y\ x > |x|oo > 9 0 . 

Here 9q is the constant that appeared in Lemma |7.51 and we used property ([6|) of 
the state space E. 

Let v > vq . We shall show that the event on the left-hand side of (1441) is contained 
in the event in (l37l) with 9 = i> 1//d . Suppose the event in (1411) happens, so that some 
x < ql satisfies 

(42) v -i/W+i) \ x \^i) > a ) _ and > v _ 

Note that the above inequality forces a(x,t) > — 00 , while the earlier requirement 
on vo forces <r(x) < 00 , and thereby also <r(x,f) < 00 . Find a maximizer y < x so 
that 

cr(x,t) = a(y) + H((y, 0), (x, t)). 

Regarding the location of y, we have two cases two consider. 

Case 1. y G [x — 2a|x|ool, x]. Let y' = x — 2a|x|ool. Then |x — y'\oo > «|?/|oo by 
virtue of a G (0,1/2). Also y' < x so the choices made above imply |?/|oo > |x|oo > 
v 1/d . 

H((y',0), (x,t)) > H((y, 0), (x, t)) = a(x,t ) - a(y) 

> I(K, a) - v~ 1/{d+1) \x\ d J^ d+1) - cr(x) 

> (/?! - 2)|x|^ (d+1) > /3(2a) d/(d+1) \x\^ (d+1) 

= /3\x-y'\ d oi {d+1) . 
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In addition to (l42l) . we used — v 1 C d + 1 ) > —1, —a(y) > — a(x ) > Pi\x\oU' d+1 \ and 
I{K,a) > -vl /{d+1) > -|®| d d (d+1) . 

Case 2. y [x — 2a\x\ x l, x\. This implies \x — y|oo > a|y|oo- 

H((y,0), (x,t)) = cr(x,t ) - a(y) > I(K,a) - v~ 1/{ - d+1) \x\ d J 0 {d+1) - a(y) 

> -\x\^ {d+1) - v~ 1/(d+1) \x\^ d+1) + Pi\y\ d J, (d+1) 

>2f3\y\ d £ d+1) >P\ x -y\ d £ d+1) - 


We conclude that the event in (14T1) lies inside the event in (1371) with 6 = v x ! d , as 
long as v > vq, and the inequality in (14T1) follows from (1371) . □ 

Corollary 7.8. Let K be a compact cube, e > 0, and 0 < t < oo. Then there exists 
a deterministic compact cube L such that 

P a {S k (<t( s )) C L for all s £ [0,f]} > 1 — e. 


Proof. For 0 < s < t, x £ Sk(ct(s )) implies that I(K,a(s)) is finite, x < ql and 
<r(x,s) > I(K,a(s)). Consequently 

l^loo — ^ q,I(K,a(s)) (°’( s )) — ■ 

Thus given e, we can choose L = [—ml, ml] with m picked by Lemma 17.71 so that 

P a {'k q j(K,*)(v(t)) > m d } < e. □ 

We are ready for the last stage of the proof of Theorem 13.11 

Proposition 7.9. Let ft be a bounded measurable function on E supported on the 
compact cube K = [— ql,ql] of R d , and a £ E. Then 

(43) E a [(f>(cr(t))\ — <f(a) = f E a [C(l)(a(s))]ds. 

Jo 

Proof. Pick a small r > 0 so that t = mr for an integer m, and denote the partition 
by Sj = jr. By the Markov property, 


E^<f(am - cf(a) = E^ 

771—1 

^{^)[0(a(r))]-0(a( S ,-))} 

_j=0 

= E a 

m-1 

/ Yl 1 (s j ,s j+1 ](s)£0(cr(s j+1 ))ds 

Jo J=0 


(44) 


+ - E a [Cf>{a(t.))}) + E a 


771—1 


A r( CT ( s :i)) 

l=o 


where the terms A T (cr(sj)) are as defined in Lemma 17.41 

We wish to argue that, as m —+ oo and simultaneously r —> 0, expression (1131) 
after the last equality sign converges to the right-hand side of (1131) . 

Note first that C(f{a) is determined by the restriction of er to the set Sk{v) OK. 
By Corollary 17.81 there exists a fixed compact set L such that Sk{&(s)) U K C L 
for 0 < s < t with probability at least 1 — e. By Corollary 17.61 the time evolution 
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{cr(x, s ) : x £ L, 0 < s < t} is determined by the finitely many Poisson points in the 
random compact rectangle \—M,M) d x [ 0 ,t]. Consequently the process C<f(a(s)) is 
piecewise constant in time, and then the integrand l(sj,sj+i]( s )£^( <T ( s .j'+i)) 

converges to Ccf>(a(s )) pointwise as m —> oo. This happens on an event with prob¬ 
ability at least 1 — e, hence almost surely after letting e —> 0. 

To extend the convergence to the expectation and to handle the error terms, we 
show that 


(45) 


E a 


sup V’ic(er(s)) 

0 <s<t 


< OO. 


Before proving (l45|) . let us see why it is sufficient. Since 


(46) \Cfi{<j)\ < 2||</>|| 00 |S'_fs-(cr)|, 

(l35l) and (l45l) imply that also the first expectation after the equality sign in (1441) 
converges, by dominated convergence. The second and third terms of (14411 vanish, 
through a combination of Lemma 17.41 (l35)l . and (l45l) . 

By the bound in (THil) for sup 0<s<t V’i:( cr ( s )) and by Lemma 17771 it only remains 
to show that 

E a [ ( J(K , a{t)) - I(K, ct( 0)) + 1) 2 ] < oo. 

This follows from property ([6]) of a and the bounds for H in Lemmas 17.11 and 17.51 
We omit the proof since it is not different in spirit than the estimates we already 
developed. □ 

This completes the proof of Theorem 13.11 


8. Proof of the limit for the height function 

Introduce the scaling into the variational formula 0 and write it as 

(47) a n (nx, nt) = sup {a n (ny, 0) + H((m/, 0), (nx, nt))}. 

y(zt{. d :y<.x 

Lemma 8.1. Assume the processes a n satisfy m and 0. Fix a finite T > 0 and 
a point b £ R d such that b > 0, and consider the bounded rectangle [— b,b\ C R ci . 
Then with probability 1 there exist a random N < oo and a random point a £ R d 
such that 

(48) & n (nx, nt) = sup {a n (ny, 0) + H((m/, 0), (nx, nt))} 

y£[a,x] 

for x £ [—6, b], t £ (0, T], n > N. 

Proof. For (3 > e 2 T 1 /( d+1 ) and b £ R d fixed, one can deduce from Lemma O and 
Borel-Cantelli that, almost surely, for large enough n, 

H((to, 0), ( nb , nt.)) < (3n\b - z|^ (d+1) 

for all i £ Z d such that i < b and |i — b > 1. If y £ R d satisfies y < b and 
| y — foloo > 1, we can take i = [y\ (coordinatewise integer parts of y) and see that 

(49) H((ra/, 0), ( nb , nt)) < fin + f3n\b - y\ d ^ d+l) 
for all such y. 




Growth model 


223 


In assumption (fl3l) choose C > (3 so that —C+(2 + \b\i£ < ' d+1 ' > )/3 < uo(—b) — 1. Let 
N and M be as given by m, but increase M further to guarantee M > 1. Now take 
a £ R rf far enough below —b so that, if y < b but y > a fails, then \y\oo > M. [Since 
assumption (1131) permits a random M > 0, here we may need to choose a random 
a £ R d .] Then by (THU) , if y < b but y > a fails, then a n (ny, 0) < —Cn\y\ < ^ d+1 \ 
Now suppose x £ [—b, b], y < x, but y > a fails. Then 

cr n (ny, 0) + H((ray, 0), (rax, nt)) 

< &n(ny, 0) + H((ray, 0), (nb, nt)) 

< -Cn\y\ d ^ d+1) +pn + (3n\b - y\ d ^ d+1) 

< n ([-C + P)\y\ d oi (d+1) + P + mt {d+1) ) 

< nuo(-b) — n < a n {—nb , 0) — ra/2 

[by assumption (fl2l) . for large enough ra] 

< cr„(rax,0) — ra/2 [by monotonicity]. 

This shows that in the variational formula 63) the point y = x strictly dominates 
all y outside [a, x]. □ 

Starting with (H51) the limit IT5l) is proved (i) by partitioning [a, x] into small 
rectangles, (ii) by using monotonicity of the random variables, and the monotonicity 
and continuity of the limit, and (iii) by appealing to the assumed initial limits 63 
and to 


(50) 


L H((ray, 0), (rax, nt)) -*■ c d+i ((x - y)\t) 1/(d+1) = tg((x - y)/t) 


a.s. 


To derive the limit in © from ([3| one has to fill in a technical step because in 
(1501) the lower left corner of the rectangle (ny, rax] x (0, nt] moves as ra grows. One 
can argue around this complication in at least two different ways: (a) The Kesten- 
Hammersley lemma j28j, page 20, from subadditive theory gives a.s. convergence 
along a subsequence, and then one fills in to get the full sequence. This approach 
was used in 24]. (b) Alternatively, one can use Borel-Cantelli if summable deviation 
bounds are available. These can be obtained by combining Theorems 3 and 9 from 
Bollobas and Brightwell [6|. 


9. Proof of the defect boundary limit 

In view of the variational equation 0, let us say cr(x, t) has a maximizer y if y < x 
and cr(x, t) = <j{y, 0) + H((y, 0), (x, t)). 

Lemma 9.1. Suppose two processes a and £ are coupled through the space-time 
Poisson point process. 

(a) For a positive integer m, let D m {t) = {x : > <j{x,t) + to}. Then if 

x £ D m (t), C(x,t) cannot have a maximizer y £ D m (0) c . And if x £ D m (t) c , a(x,t) 
cannot have a maximizer y £ D m { 0). 

(b) In particular, suppose initially <j(y, 0) < f(y, 0) < a(y, 0) + h for all y £ R d , 
for a fixed positive integer h. Then this property is preserved for all time. If we 
write 

Aft) = {x : £(x, t) = <r(x, t) + h}, 

then 

(51) Aft) = (x : a(x,t) has a maximizer y £ A(0)}. 
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(c) If h = 1 in part (b), we get additionally that 
(52) Aft) c = {x : £(x,t) has a maximizer y £ A(0) c }. 


Proof, (a) Suppose x £ D m (t), y £ D m ( 0) c , and £{x,t) = £(y, 0) + H((y, 0), (x,f)). 
Then by the definition of D m (t ), 

t) <C(x,t)-m = c (y, o) - TO + H ((y, 0), (x, t)) < a(y , 0) + H ((y, 0), (a:, t)) - 1 

which contradicts the variational equation 0. Thus £(x, t ) cannot have a maximizer 
y £ D m { 0) c . The other part of (a) is proved similarly. 

(b) Monotonicity implies that cr(x,t) < £(x,t) < a(x,t) + h for all (x,f), so 
A(t) = Dh(t). Suppose x £ Aft). By (a) £(x,t) cannot have a maximizer y £ A(0) c , 
and so £{x,t) has a maximizer y £ A(0). Consequently 

x, t) = £(x, t) - h = £{y, 0) - h + H((y, 0), (x, t)) = a{y , 0) + H((y, 0), (x, t)), 

which says that <r(x, t ) has a maximizer y £ A(0). On the other hand, if <r(x, t) has 
a maximizer y £ A(0), then by (a) again x £ A(t) c . This proves (1511) . 

(c) Now Aft) = D\{t) and A{t) c = {x : a(x,t) = f(x,t)}. If (f(x,t) has a 

maximizer y £ A(0) c , then by part (a) x ^ Aft). While if x £ A(t) c , again by part 
(a) er (x,t) must have a maximizer y £ A(0) c , which then also is a maximizer for 
C(x,t). This proves (l52l) . □ 

Assume the sequence of processes a n (■) satisfies the hypotheses of the hydrody¬ 
namic limit Theorem 14.11 which we proved in Section [8] The defect set A n ft) was 
defined through the (a n ,£ n ) coupling by (l27l) . By (HTT1) above, we can equivalently 
define it by 

(53) A n (t) = {x : <j n (x,t) has a maximizer y £ A„(0) }. 

In the next lemma we take the point of view that some sequence of sets that depend 
on at has been defined by (1531) . and ignore the (<r n , £„) coupling definition. 

Lemma 9.2. Let B C R d be a closed set. Suppose that for almost every sample 
point uj in the underlying probability space, a sequence of sets A n ( 0) = A n (0;o;) is 
defined, and has this property: for every compact K C R tZ and e > 0, 

(54) {n -1 A„(0)} (~l K C (~l K for all large enough n. 

Suppose the sets A n ft) satisfy (1551) and fix t > 0. Then almost surely, for every 
compact K C R d and e > 0, 

(55) {n -1 A n (nf)} fl K C W(B, t)^ D K for all large enough n. 

In particular, ifW(B,t) = 0, then (1551) implies that {n~ 1 A n (nt)} D K = 0 for all 
large enough n. 

Proof. Fix a sample point u> such that assumption (1541) is valid, the conclusion of 
Lemma f8. II is valid for all b £ Z^, and we have the limits 


(56) 


n 1 (j n fnx,nt) —* u(x,t) for all (x,t), 
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and 


(57) n 1 H((nj/, 0), (nx, nt)) —* tg((x — y)/t) for all y, x, t. 

Almost every to satisfies these requirements, by the a.s. limits (15U1) and US]), by 
monotonicity, and by the continuity of the limiting functions. It suffices to prove 
(1551) for this fixed lo. 

To contradict (1551) . suppose there is a subsequence rij and points Xj £ K 
such that rijXj £ A nj (rijt) but Xj ^ W(B,t)( e \ Note that this also contradicts 
{n _1 A n (nt)} n K = 0 in case W(B,t) = 0, so the empty set case is also proved by 
the contradiction we derive. 

Let rijUj £ A nj (0) be a maximizer for a nj ( rijXj, rijt ). Since the Xj’s are bounded, 
so are the yj’s by Lemma 18.11 and we can pass to a subsequence (again denoted 
by {j}) such that the limits Xj — ► x and yj —► y exist. By the assumptions on Xj, 
x 0 W(B,t). For any e > 0, yj £ B^ for large enough j , so y £ B by the closedness 
of B. 

Fix points x' < x" and y' < y" so that x' < x < x" and y' < y < y" in the 
partial order of R d . Then for large enough j, x' < Xj < x" and y' < yj < y". By 
the choice of yj, 

a nj (njXj,rijt) = a nj (rijyj, 0) + (njXj, njt)) 

from which follows, by the monotonicity of the processes, 

n J 1(T n j {rijx ', rijt) < nJ 1 <j nj (rijXj, rijt) 

< n~ 1 a nj {rijy" , 0) + n~ l H{{rijy' , 0), (njx”, rijt)). 

Now let rij —> oo and use the limits (|5fil) and m to obtain 

u(x', t) < uo{y") + tg((x" - y')/t). 


We may let x',x" —> x and y',y" —> y, and then by continuity u(x,t) < uo(y) + 
tg((x — y)/t). This is incompatible with having x W(B,t) and y £ B. This 
contradiction shows that, for the fixed u>, (1551) holds. □ 

We prove statement (l29l) of Theorem 15.11 The assumption is that 
(58) B n K C |n _1 A n (0)} D K C B^ D K for all large enough n. 

We introduce an auxiliary process £ n (x,t ). Initially set 


(59) 


€n{y, 0) 


<7n{y, o), y £ A n (0) 
&n{y, 0) + 1, y £ A„(0). 


£„(?/, 0) is a well-defined random element of the state space E because A„(0) is 
defined m in terms of Cn(y, 0) which lies in E. Couple the process £„ with a n and 
C, n through the common space-time Poisson points. Then 


a n {x,t) < £ n {x,t) < cr n {x,t) + 1. 

By part (b) of Lemma liTTl A n (t) that satisfies (1551) also satisfies 
An(t) {X . ^ n (^X,t) (7 n {x,t) 1). 
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Then by part (c) of Lemma 19.11 

(60) A n (t) c = {x : £ n (x,t) has a maximizer y G A n ( 0) c }. 

The first inclusion of assumption (l58l) implies that n~ 1 A n (0) c r\K C (B c )^ C\K 
for large n. The processes inherit all the hydrodynamic properties of the processes 
a n . Thus by (|fiUl) we may apply Lemma T9.21 to the sets A n (nt) c and the processes 
£ n (nt) to get 

(61) n~ 1 A n (nt) c n K C W(W,t ) (a) n K 
for large enough n. By (l55jl and (RTTTl . 

bd {n- 1 A n (nt)}r\K C W{B, t)W n W(W, t) {5) n K 

for large n. For small enough S > 0, the set on the right is contained in 
[W(B,t) nW(T^,t)] {E) HK = X(B,i)& nK. This proves ([29]l . 


To complete the proof of Theorem 15.11 it remains to prove 
(62) 

W (B, t)^ ^ D K C n 1 A n (nt ,) (~l K C W(B , t)^ fl K for all large enough n 

under the further assumption that no point of W(B c ,t) is an interior point of 
W(B,t). 

The second inclusion of (RSI) we already obtained in Lemma SSI (RJT1) implies 


W( B c , t) (,5) l C n K C nr'Anint) n K. 


It remains to check that, given e > 0, 


kF(B,t) ( " e) 


n K c 


W(B C , t) 


(<5) 


nx 


for sufficiently small S > 0. Suppose not, so that for a sequence Sj \ 0 there exist 
Xj G W(B,t)^ e ^ fl W( B c , t)^ 5 A D K. By Lemma [6.11 the set W(B c ,t) is closed. 
Hence passing to a convergent subsequence Xj —> x gives a point x G W(B c ,t) 
which is an interior point of W(B,t), contrary to the hypothesis. 


10. Technical appendix: the state space of the process 

We develop the state space in two steps: first describe the multidimensional Skoro- 
hod type metric we need, and then amend the metric to provide control over the 
left tail of the height function. This Skorohod type space has been used earlier (see 
[5) and their references). 

10.1. A Skorohod type space in multiple dimensions 

Let (X,r) be a complete, separable metric space, with metric r(x,y) < 1. Let 
D = D(R d , X ) denote the space of functions o : R d —► X with this property: for 
every bounded rectangle [a, b) C R d and e > 0, there exist finite partitions 

at = s° < s] < ■ ■ ■ < sT = bi 
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of each coordinate axis (1 < i < d) such that the variation of cr in the partition 
rectangles is at most e: for each k = (fci, k 2 , ■ ■ ■ , fed) £ nf=i{ 0 , 1 , 2 ,..., m* — 1 }, 

(63) sup{r(cr(a;), cr(y)) : < Xi,yi < s ^ i+1 (1 < i < d)} < e. 

Note that the partition rectangles are closed on the left. This implies that cr is 
continuous from above: a(y ) —> cr(x) as y —> x in R d so that y > x; and limits exist 
from strictly below: lim cr{y) exists as y —> x in R d so that y < x (strict inequality 
for each coordinate). 

We shall employ this notation for truncation in R d : for real u > 0 and x = 

(xi,...,x d ) £ R d , 

[x\ u = ({xi A«)V (-it), {x 2 Au)V (~u ),..., (xd A u) V (-u)). 

Let A be the collection of bijective, strictly increasing Lipschitz functions A : R d —> 
R d that satisfy these requirements: A is of the type A(aq,... ,Xd) = (Ai(xi),..., 
A d(xd)) where each A,; : R —* R is bijective, strictly increasing and Lipschitz; and 

7U) = 7o (A) +7i(A) < oo 

where the quantities 70 (A) and 71 (A) are defined by 


d 

7 o(A) = V sup 

s,££R 


log 


A j(t) - Xj{s) 
t — s 


a nd 

71 (A) = / e -u (l A sup I [A(x)]„ - [x\ u I ) du. 
Jo v x€R d J 

For p, cr S D, A € A and u > 0, define 


d{p,cr, X,u) = sup r[p(\x] u ) , cr([A(x)] tl )). 
xeR d 


And then 
(64) 


^s(p,cr) = inf 
AeA 


7 (A) V / e u d(p , cr, A, u) du 


The definition was arranged so that 7 (A -1 ) = 7 (A) and 7 (A o p) < 7 (A) + 'y(p), so 
the proof in llj, Section 3.5, can be repeated to show that ds is a metric. 

It is clear that if a sequence of functions cr„ from D converges to an arbitrary 
function cr : R d —> A', and this convergence happens uniformly on compact subsets 
of R d , then a £ D. Furthermore, we also get convergence in the c?s-metric, as the 
next lemma indicates. This lemma is needed in the proof that (D,ds) is complete. 


Lemma 10.1. Suppose a n , a € D. Then ds(cr„, a) 0 iff there exist A„ g A such 
that y(A„) —> 0 and 

r(a n (x),a{X n (x))) —> 0 
uniformly over x in compact subsets of R d . 


Proof. We prove ds(cr„, cr) —> 0 assuming the second condition, and leave the other 
direction to the reader. For each rectangle [—Ml, Ml), M = 1, 2,3,... , and each 
e = 1/A', K = 1, 2, 3,... , fix the partitions {s/} that appear in the definition (Ri51) 
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of a £ D. Pick a real u > 0 so that neither u nor —u is among these countably 
many partition points. 

d(a n ,a, X n , u) = sup r(a n ([x ] u ), <j([A„(a;)]„)) 

iceR d 

< sup r(cr n ([x ] u ), cr(X n ([x] u ))) + sup r(cr(A„([x] u )), cr([A„(a:)] u )). 

The first term after the inequality vanishes as n —> oo, by assumption. 

Let e = 1/K > 0, pick a large rectangle [—Ml, Ml) that contains [— ul,ul] 
well inside its interior, and for this rectangle and e pick the finite partitions that 
satisfy J63| for a, and do not contain ±tt. Let d > 0 be such that none of these 
finitely many partition points lie in (±n — S, ±it + <5). If n is large enough, then 
sup xe r_ M1 ]Afl ] A„(;r) — x\ <5, and one can check that A„([a:] u ) and [A„(a;)] w lie in 
the same partition rectangle, for each x £ R d . Thus 


sup r (cr(A„([a;] u )), cr([A n (a;)] u )) < e. 

x£R d 


We have shown that d(a n , a, X n ,u) —> 0 for a.e. u > 0. 


□ 


With this lemma, one can follow the proof in [11], page 121, to show that ( D , ds ) 
is complete. Separability of ( D , ds) would also be easy to prove. Next, we take this 
Skorohod type space as starting point, and define the state space E for the height 
process. 


10.2. The state space for the height process 

In the setting of the previous subsection, take S = Z* = ZU{±oo} with the discrete 
metric r( x, y) = l{a: fy y}. Let E be the space of functions a £ D( R d , Z*) that are 
nondecreasing [a(x) < cr(y) if x < y in R d ] and decay to —oo sufficiently fast at 
—oo, namely 

(65) for every b £ R d , lim sup { \y\^^ d+1 ^a(y) : y <b, \y\oo > m\ = -oo. 

Condition (1651) is not preserved by convergence in the ds metric, so we need to fix 
the metric. 

For fj £ E, h £ Z, and b £ R d , let y b,h (a) be the maximal y < b in R d such that 
the rectangle [y, b] contains the set {x <b : a(x) > h}. Condition (l65l) guarantees 
that such a finite y b,h (a) exists. In fact, (R)5ll is equivalent to 

(66) for every b £ R d , lim \h[’< d+r > /d \y b ’ h {a)\ 00 = 0. 

h —>—oo 

For p, a £ E and b £ R d , define 

9 b {p,<j) = sup |/ify (d+1)/d • | y b ' h {p) -y b ' h (cr )|oo 

h<-l 

and 

0(p, cr) = f e~ |b|oo (lA 6 b (p,a))db. 

J R d 

0(p, a) satisfies the triangle inequality, is symmetric, and 0(cr, a) = 0, so we can 
define a metric on E by 

ds(p, cr) = 0(p, a) + d s (p , cr). 

The effect of the 0(p, a) term in the metric is the following. 
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Lemma 10.2. Suppose ds{<J n , cr) —> 0. Then d E (cr„, cr) —> 0 iff for every h £ R d , 
(67) lim sup |/ l |- ( rf+ 1 ) / rf|^^( CT „)| oo = o, 

h —>• —oo n 


or equivalently, for every h £ R d 


( 68 ) 


lim sup 

M—> oo n 


sup 

y<b 

I 9|=o>M 


°n(y) 

\y& d+1) 


—oo. 


We leave the proof of the above lemma to the reader. Lemmas 110.11 and 110.21 
together give a natural characterization of convergence in E. 

Lemma 10.3. The Borel a-field By: is the same as the a-field T generated by the 
coordinate projections a i—> <j(x). 

Proof. The sets {x : a(x) > h} are closed, so the functions a i—> a(x) are upper 
semicontinuous. This implies T C By- 

For the other direction one shows that for a fixed p £ E, the function a i—> 
d E (p,cr) is jF-measurable. This implies that the balls {a £ E : d E (p, cr) < r} are 
^-measurable. Once we argue below that E is separable, this suffices for By Q F■ 
To show the ^-"-measurability of a i—> ds ( p , cr) one can adapt the argument from 
page 128 of 11]. To show the ^-measurability of cr i—> 0(p, cr), one can start by 

y b ’ h (a) from R d x E 

□ 


arguing the joint £>Rd ® iF-measurability of the map (6, cr) 
into R d . We leave the details. 


The remaining work is to check that (E, d E ) is a complete separable metric space. 

Proposition 10.4. The space (E, d E ) is complete. 

We prove this proposition in several stages. Let {cr„} be a Cauchy sequence in 
the d E metric. By the completeness of ( D,ds ), we already know there exists a 
a £ £>(R d ,Z* ) such that ds(cr„,cr) —> 0. We need to show that (i) a £ E and (ii) 
0(cr„,cr) —> 0. 


Following the completeness proof for Skorohod space in [11], page 121, we may 
extract a subsequence, denoted again by cr„, together with a sequence of Lipschitz 
functions ip n £ A (actually labeled in ll]|), such that 


(69) 

and 


7 W>n) < 2 1 


(70) 


cr n {ipn{x)) —■” cr(a;) uniformly on compact sets. 


Step 1. cr £ E. 

Fix & £ R d , for which we shall show (Rhjl) . It suffices to consider b > 0. Let b k = 
b + k 1. By passing to a further subsequence we may assume 0(cr„, a n +i) < e~ n . 
Fix no so that 

(71) exp(|& 2 |oo + d{n + 1) — n 2 ) < 2~ n for all n > no- 

Lemma 10.5. For n > no there exist points f3' n in R d such that b 1 < (3 n+1 < f3 n < 
b 2 , and 0pn(a n ,cr n+1 ) < 2~ n . 
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Proof. Let a n = b 1 + e n ■ 1 in R d . 

2 

6 P 0((J n , ) 

> inf {l/\9 x (cr n ,a n+ i)} ■ e~\ b \°° ■ Leb d {x : a n+ i < x < a„}, 

xe(a n + l,OLn) 

where 

Lebd{a; : a n+ \ < x < a n } = ( e~ n - e~ n ~ 1 ) d > e _d («+i) 

is the d-dimensional Lebesgue measure of the open rectangle (a n +i,a n ). This im¬ 
plies it is possible to choose a point (3 n G (a n +i,a n ) so that (oVi, o'™ +i) < 

2~ n . □ 

/3 n+1 < /3 n implies j/ 3 ” +1 > h (<j n +i) > yP n ' h (a n +i) — (/ 3 n — fi n+1 ). For each hxed 
h < — 1, applying the above Lemma inductively gives for n > rig: 

yP n+1 ’ h (a n+1 ) > y? n ’ h (cT n+1 ) - (/?” - (3 n+1 ) 

> y pn ' h {(Tn) ~ |ft| (d+1)/d 2- n • 1 - (P n - f3 n+1 ) 

n 

>■■■> /"’■>„„) - \h\( d+1 ^ d 2 ~ k ■ i - (Z 3 ” 0 - z 3 ”* 1 ), 

k=no 

from which then 

(72) inf y bl ’ h (a n ) > - | h \(d+D/d 2 i-n 0 . 1 _ ^ ^ 

n>n o 

Now fix h < —1 for the moment. By m we may fix a rectangle [y 1 ,^ 1 ] that 
contains the sets {x < b l : a n (x) > h} for all n > uq. Let Q = [y 1 — l,b x + 1] 
be a larger rectangle such that each point in [j/ 1 , b 1 ] is at least distance 1 from 
Q c . By (IPPll and (171)1) we may pick n large enough so that \tp n {x) — x\ < 1/4 and 
o’n{'4’n{x)) = &(x) for x G Q. [Equality because Z* has the discrete metric.] 

We can now argue that if x < b and a(x) > h , then necessarily x G Q, ip n (x) < b 1 
and a n (tp n (x)) > h , which implies by d72l) that 

x > if n (x) - (1/4)1 > y bl ’ h (a n ) - (1/4)1 > /’^Ko) - (5/4 + \ h \id+D/d 2 i-n 0) . x 
This can be repeated for each h < —1, with no fixed. Thus for all h < —1, 

\y bd \a)\ < |&| V (\y b2 ’ h ((r no )\ +5/4+ |/i|( d+1)/d 2 1 - n o) , 
and then, since a no G £, 

lim \h\~( d+1) / d \y b ’ h (a)\ < 2 1 ~ n °. 

h—> — oo 

Since no can be taken arbitrarily large, (l66l) follows for <r, and thereby a G £. 
Step 2. Q(a n ,a) —> 0. 

As for Step 1, let us assume that we have picked a subsequence a n that satisfies 
(1001) and (1701) and 0(cr n+ i, a n ) < e - ™ 2 . Let />„ = , 0“ 1 . If we prove 0(er n ,cr) —> 0 
along this subsequence, then the Cauchy assumption and triangle inequality give it 
for the full sequence. 
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Fix an arbitrary index rq and a small 0 < £o < 1. Fix also [3 g R d . For each 
h < —1, fix a rectangle [y h ,/3] that contains the sets {x < (3 : cr n (x) > h} for each 
a n for n > n i, and also for cr, which Step 1 just showed lies in E. This can be done 
for each fixed h because by d72l) there exists no = no{(3) defined by (1711) so that 
the points are bounded below for n > no- Then if necessary decrease y h 

further so that 


y h < /^( CTl ) a y 0 ' h (a 2 ) A • • • A y P ' h {a no - 1 ). 

Let Qh,k = [ y h — k 1, (3 + fcl] be larger rectangles. 

On the rectangles Qh : 2 , h < —1, construct the finite partition for cr which sat¬ 
isfies (Rxil) for e = 1/2, so that the discrete metric forces cr to be constant on the 
partition rectangles. Consider a point b = (bi, b 2l ..., bd) < (3 with the property 
that no coordinate of b equals any one of the (countably many) partition points. 
This restriction excludes only a Lebesgue null set of points b. 

Find S\ = E\(j3,b,h) > 0 such that the intervals (bi — e±,bi +£i) contain none 
of the finitely many partition points that pertain to the rectangle Qh, 2 - Pick n = 
n(/3,b,h) > n\ such that cr„ (ip n (a:)) = °’(x) and | ip n (x) — x\ < (&o A £i)/ 4 for 
x G Qh,2- Since the maps ?/>,</> do not carry any points of [y h ,f3] out of Qh,i , 
y h ’ b {a n ) = y h,b (a o <j> n ). It follows that 

I/’V) - /’Vn)I = I/’V) - /’V O K)\ < so. 

The last inequality above is justified as follows: The only way it could fail is that cr 
(or cr o <f> n ) has a point x < b with height > h, and cr o </>„ (respectively, cr) does not. 
These cannot happen because the maps i[>, <j> cannot carry a partition point from 
one side of hi to the other side, along any coordinate direction i. 

Now we have for a.e. b < (3 and each h < — 1, with n = n{(3 , b , h) > ny. 

\ h \-(d+i)/d\ y h ,b( a j _ y h ’ b (a ni )\ 

< \h\- ( - d+1 V d \y h ’ b (a) - y h ’ b (a n )\ + \h\^ d+1 ^ d \y h ’ b (a n ) - y h ’ b (a ni )\ 

— “1“ Qb^P’ni ^ni) 

< £o + Sup 0 b (a m i &n\ ) • 

m:m>ni 

The last line has no more dependence on (3 or h. Since (3 was arbitrary, this holds 
for a.e. b G R d . Take supremum over h < — 1 on the left, to get 


1 A 0 b (cr, a ni ) < e 0 + sup {1 A 9 b {a m , a ni )} for a.e. b. 

m:m>n i 


Integrate to get 


Jn d 


<Ce 0 + [ e~^°° 

sup (1 A 0f,(cr m , cr ni )} 

J R d 


<Ce o + [ e- |b| “ 

f m—1 

sup < ^ 1 A0 b (cr fc+ i,c 

Jn d 

^ k—ni 

oo „ 

= Ce o + / e 

fc=m jRd 

{l /\ 9 b (ak+i,(Xk)} db 

OO OO 

= Ceq + ©(crfc_|_i,crfc) < C^o + e , 


k=n± 


k=ni 
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where C = f Rd e db. Since n i was an arbitrary index, we have 

lim sup 0 (cr ni , a) < Ceo- 

n± —kx> 

Since eo was arbitrary, Step 2 is completed, and Proposition 110.41 thereby proved. 

We outline how to construct a countable dense set in (E ,ds). Fix a < b in 
Z d . In the rectangle [a, b] C R d , consider the (countably many) finite rational 
partitions of each coordinate axis. For each such partition of [a, b] into rectangles, 
consider all the nondecreasing assignments of values from Z* to the rectangles. 
Extend the functions a thus defined to all of R rf in some fashion, but so that they 
are nondecreasing and Z*-valued. Repeat this for all rectangles [a, b\ with integer 
corners. This gives a countable set D of elements of D(R d , Z*). Finally, each such 
a G D yields countably many elements a G E by setting 

Z( T ) = I” 00 ’ < h 

( cr(x), a(x) > h 

for all h £ Z. All these a together form a countable set ECE. 

Now given an arbitrary a £ E, it can be approximated by an element a £ E 
arbitrarily closely (in the sense that a = a o (f> for a map <j> £ A close to the identity) 
on any given finite rectangle [—/?,/?], and so that y b,h (a) is close to y b,h (a) for all 
b in this rectangle, for any given range ho < h < —1. Since < £ 

for h < ho for an appropriately chosen ho, this suffices to make both d{a,a,4>,u) 
and 9b(cr,d) small for a range of u > 0 and b e R d . To get close under the metric 
ds it ^suffices to approximate in a bounded set of u' s and V s, so it can be checked 
that E is dense in E. 
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